Rate of equilibration of a one-dimensional Wigner crystal 
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Abstract 

We consider a system of one-dimensional spinless particles interacting via long-range repulsion. In the limit of strong 
interactions the system is a Wigner crystal, with excitations analogous to phonons in solids. In a harmonic crystal 
the phonons do not interact, and the system never reaches thermal equilibrium. We account for the anharmonism of 
the Wigner crystal and find the rate at which it approaches equilibrium. The full equilibration of the system requires 
umklapp scattering of phonons, resulting in exponential suppression of the equilibration rate at low temperatures. 
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The low-temperature physics of interacting electron 
systems is usually described in the framework of the so- 
called Luttinger liquid theory [1] . The phenomenologi- 
cal nature of this approach enables one to study the sys- 
tems with any interaction strength, provided that the 
physics is controlled by the low-energy excitations. On 
the other hand, even at low temperature T some phe- 
nomena involve excitations with energies much higher 
than T. In such cases microscopic approaches are usu- 
ally more effective. 

An example of such a phenomenon is the equilibra- 
tion of a one-dimensional system. The latter is under- 
stood most easily in the case of weak interactions be- 
tween electrons, when the usual picture of quasiparticle 
excitations is applicable. Because of the conservation 
of momentum and energy in electron-electron scatter- 
ing, the two-particle collisions do not affect the distri- 
bution function. Thus at weak interactions the equi- 
libration is accomplished via the three-particle colli- 
sions, Fig. 1(a). To equilibrate the chemical potentials 
of the right- and left-moving electrons, the collisions 
must change the numbers of electrons on each branch, 
i.e., the backscattering events, such as the one shown in 




Corresponding author. E-mail: matveev@mailaps.org 



2k F q 



Fig. 1. (a) Three particle scattering process leading to equi- 
libration of weakly-interacting one-dimensional electrons, (b) 
Spectrum of the hole excitation. Backscattering of the hole 
occurs when its wave vector crosses the point q — kp. 

Fig. 1(a), are needed. The most effective such process 
includes backscattering of a particle at the very bot- 
tom of the band [2]. The small probability of finding 
an available final state deep below the Fermi level Ef 
results in exponential suppression of the equilibration 
rate t" 1 oc e~ EF/T [3]. 

The key event in the process of equilibration is the 
backscattering of a hole at the bottom of the band. It 
is therefore helpful to focus on the motion of such a 
hole in momentum space. For electrons with quadratic 
spectrum Ek = fi 2 k 2 /2m the energy of the hole exci- 
tation is e{q) — Ek F — Ek F - q = hvirq(l — q/2ki?), as- 
suming that the hole is created by moving an electron 
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from state JtF — q to the right Fermi point fcf- Here Uf 
is the Fermi wave vector and vp = hkp/m is the Fermi 
velocity. Backscattering occurs when the wave vector 
of the hole q crosses Jcf, Fig. 1(b). 

Because the equilibration of the system involves ex- 
citations with energies of the order of the bandwidth, 
this phenomenon is not captured by the conventional 
Luttinger liquid theory [1]. This makes the generaliza- 
tion of the above picture beyond the weakly- interacting 
limit rather challenging. On the other hand, it is possi- 
ble to develop a microscopic theory of this phenomenon 
in the case of strong Coulomb interactions [4] . Here we 
revisit this approach and obtain the full expression for 
the equilibration rate for arbitrary strong long-range 
repulsion. The key idea is that as repulsion of electrons 
becomes stronger, the system minimizes its energy by 
forming a periodic structure known as the Wigner crys- 
tal. Although the long-range order in such a system 
is destroyed by quantum fluctuations [5], the presence 
of the strong short-range order enables us to treat the 
system as anharmonic chain described by the Hamil- 
tonian 
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Here pi and xi are the momentum and coordinate of 
the l-ih particle and V(x) is the interaction potential. 

The excitations of the system are essentially the 
phonons in the electronic crystal. They are conve- 
niently described in terms of the displacements ui = 
xi — la of electrons from their equilibrium positions, 
where a is the mean interparticle distance. Strong 
repulsion means small displacements, \ui — -C 
\l — l'\a, and the Hamiltonian can be approximated by 
that of a harmonic chain 
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where the r-th derivative of V(x) is denoted as 
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The phonon modes of this Hamiltonian are easily 
found, 
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At small wave vector qa — > the excitation spectrum 
is linear, w q — s\q\, where s = I 2 /m) 1//2 is the 

"sound velocity" in the Wigner crystal measured in 
units of lattice spacings per unit time. The spectrum 
is periodic in q, with the Brillouin zone — it < q < n, 
Fig. 2 




Fig. 2. Excitation spectrum of a onc-dimcnsional Wigner crys- 
tal. Equilibration processes involve a phonon qi being scat- 
tered to a state q 1 outside the Brillouin zone as a result of 
a collision with a thermal phonon qi . Conservation of energy 
and momentum requires the latter to be backscattcrcd to q'^. 
This umklapp process is analogous to a hole crossing the point 
q — kp in the weakly interacting Fermi gas, Fig. 1(b). (Note 
that if the distances arc measured in units of a, the Fermi wave 
vector kp — 7T.) 

Any eigenstate of the harmonic chain (2) can be de- 
scribed by the set of occupation numbers iVq of the 
phonons. In the absence of interaction of phonons, the 
lifetime of any such state is infinite, and the system 
never reaches thermal equilibrium. On the other hand, 
the harmonic Hamiltonian (2) is merely the leading 
term of the expansion of Eq. (1) in the small param- 
eter K — nh/ma 2 s. The next term is proportional to 
the third power of the displacement ui and generates 
scattering processes involving three phonons. At low 
temperature T <^ hw-^ the typical quasimomenta of 
phonons are small, q <g T/hs, umklapp scattering is 
suppressed, and the phonons remain in the first Bril- 
louin zone upon scattering. This means that apart from 
energy, collisions conserve the total quasimomentum of 
the phonons. This yields the equilibrium phonon dis- 
tribution 
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characterized by two parameters, the temperature T 
and the velocity u of the phonon gas with respect to 
the crystal. 

On the other hand, even at low temperatures there 
are rare umklapp collisions of phonons, such as the one 
shown in Fig. 2, which do not conserve their quasimo- 
mentum. As a result, one expects the velocity u to re- 
lax gradually as u — —u/t with a small relaxation rate 
t -1 . To find it we notice that the most efficient umk- 
lapp processes involve a phonon qi near the boundary 
q — tv of the Brillouin zone colliding with a thermal 
phonon with q2 ~ T/hs, Fig. 2. At low temperature the 
resulting change of quasimomentum \q[ —qi\ ~ T/hs is 
small compared to the typical scale qr ~ (T/h\w" I) 1 / 2 
of the distribution (5) near the edge of the Brillouin 
zone. (Here uig = d 2 uj q .) Thus the high-energy phonon 
performs a slow diffusive motion in the momentum 
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space, and its distribution function N q (t) obeys the 
Fokker-Planck equation 
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has the meaning of the diffusion constant in momentum 
space and W q , q+ s q is the rate at which a phonon q 
changes its wave vector by 5q as a result of collisions 
with other phonons. 

The Fokker-Planck equation should be solved with 
the boundary conditions 
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obtained by extending the distribution (5) beyond the 
first Brillouin zone. Such solution [4] gives the relax- 
ation law u = — u/t with the rate 
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where B = B(tt). The temperature dependence of the 
relaxation rate is dominated by the exponentially small 
probability of the occupation of phonon states near 
the edge q — n of the Brillouin zone. Expression (9) 
is analogous to the result r _1 oc e~ Ep / T for weakly- 
interacting electrons. The strong interactions between 
electrons renormalize the activation temperature from 
Ef to tko-n- in Eq. (9). 

The temperature dependence of the prefactor in 
Eq. (9) is determined by that of the diffusion constant 
B and by T~ 7 ^ 2 explicitly present in (9). The former 
can be deduced phenomenologically [6] by treating the 
phonon near q — n as a mobile impurity in a Luttinger 
liquid, for which B is known [7] to scale as 
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We therefore conclude that the equilibration rate scales 
with temperature as r" 1 oc T^ 2 e~ Hwir/T . 

The constant \ m Eq. (10) has to be determined by 
microscopic evaluation of the scattering rate W q , q +s q 
in Eq. (7). The dominant scattering process, illustrated 
in Fig. 2, involves two phonons in both the initial and 
final states. Such scattering can be accomplished ei- 
ther in the first order in four-phonon scattering ampli- 
tude or in the second order in three-phonon scattering 
amplitude. The resulting expression for the scattering 
rate has the form [4] 
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Here N is the total number of particles in the system 
and 
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with the functions fa and defined as 
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- cos(q 3 l) - cos[(qi + g 2 + q3)l] 
+ cos[(gi + q 2 )l] + cos[(qi + q 3 )l] 
+ cos[(g 2 + q-A)l}}- 

The expressions (11) and (12) are valid for any g2 
and q' 2 . At low temperature their magnitudes are small, 
M, |g 2 | < T/hs. Thus to find \ in Eq. (10) we need 
to expand A in powers of g2 and q 2 . Carrying out such 
expansion and taking into account conservation of mo- 
mentum and energy we find 
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where 8q — q 2 — q' 2 is the small momentum change as 
a result of scattering, n = 1/a is the particle density, 
v — s/n and v q — uj' q /n are the physical velocities of 
the phonons with wave vectors and q, respectively. 
To leading order in 5q we have replaced qi,q[ — > q. 

Equations (11), (13), and (14) determine the scat- 
tering rate W q , q+ s q for any q. It is easy to see that 
W q , q+ s q oc (<5g) 2 , and Eq. (7) immediately gives the 
temperature dependence (10). To find the coefficient x 
we set q = tt and obtain 
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Remarkably, \ is fully determined by the phonon spec- 
trum and its dependence on the particle density n. 

Equations (9), (10), and (15) give the complete ex- 
pression for the equilibration rate of a one-dimensional 
Wigner crystal. In the case of pure Coulomb repulsion 
V(x) — e /\x\ the velocity s diverges, and our treat- 
ment is inapplicable. However, in the experimental re- 
alizations of one-dimensional Wigner crystal, there is 
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usually a metal gate screening the interactions at large 
distances. In this case we obtain 
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where r\ — 637r 3 £(3)\/hi 2/80y / 27r is a numerical pref- 
actor, as = ft 2 /me 2 is the Bohr's radius, d is the dis- 
tance to the gate, and 
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A non-trivial test of our result (15) can be performed 
by considering the interaction potential 
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It is well known [8] that the model (18) is integrable, 
i.e., it has an infinite number of integrals of motion. 
As a result the excitations of the system have infinite 
lifetimes, and one expects the diffusion constant B to 
vanish. Our approach applies only to the limit of strong 
repulsion 7 — > 00, when the Wigner crystal approxima- 
tion is applicable. On the other hand, the parameter 
c can take any value. It is easy to obtain analytic ex- 
pressions for the phonon spectrum in the limiting cases 
c <§C n and c»n. 

At c « n one can approximate (18) with V(x) — 
y/(cx) 2 , the so-called Calogero-Sutherland model [8]. 
Then from Eq. (4) one finds 



energy given by the Debye frequency uj^ of the Wigner 
crystal. The prefactor can be expressed in terms of 
the phonon spectrum using Eqs. (9), (10), and (15). 
In the case of Coulomb repulsion the result is given by 
Eqs. (16) and (17). Finally, we have checked that the 
equilibration rate vanishes for the integrable model of 
particles with interactions in the form (18). 
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and v is given by v q at q — 0. Substitution of Eqs. (19) 
and (20) into (14) gives T q = 0, and therefore B = 0. 

At c > n the interactions fall off very rapidly with 
the distance, V(x) = ^e~ 2cx . In this case only the 
interaction of the nearest neighbor particles in the 
Wigner crystal needs to be taken into account (Toda 
lattice). The spectrum takes the form 
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As expected, substitution of Eqs. (21) and (22) into 
(14) gives T q = 0. Finally, we have checked numerically 
that the expression (15) vanishes for potential (18) for 
any c. 

To summarize, we have obtained the equilibration 
rate t -1 of one-dimensional system of particles with 
strong long-rage repulsion. At low temperatures the 
rate is exponentially suppressed with the activation 
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